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' Abstract. We investigate ciiarged and neutral current neutrino-induced incoiierent pion production off nuclei within the 

CN ■ GiBUU model at energies relevant for the MiniBooNE and K2K experiments. Special attention is paid to the entanglement 

, between measured CCQE and CCl;r+ cross sections. We further give predictions and compare to recent data measured at 

MiniBooNE. 
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INTRODUCTION 

Q . A proper understanding of neutrino induced pion production is essential for the interpretation of cuiTent neutrino 

^ ' oscillations experiments. A good knowledge of the neutrino energy is required for a precise determination of oscillation 

P5 , parameters in v^i disappearance measurements. These experiments search for a distortion in the neutrino flux at the 

detector positioned far away from the source. By comparing both, near and far neutrino energy spectra, one gains 
information about the oscillation probability and with that about mixing angles and mass squared differences. However, 
^ ' the neutrino energy is not measurable directly but has to be reconstructed from the reaction products. Present oscillation 

OO \ experiments use the CCQE reaction both as signal event and to reconstruct the neutrino energy from the outgoing muon 

CO . with two-body kinematics assuming the target nucleon is at rest. CCQE is defined as v^n — * l^p on a single nucleon; 

in the nucleus, CCQE is masked by final-state interactions (FSI). Thus, the correct identification of CCQE events is 
directly related to the question of how FSI influence the event selection. The main background to CCQE is CCl7r+ 
0^ ' production. If the pion is absorbed in the nucleus and/or not seen in the detector, these events can be misidentified as 

, CCQE. Consequently, a proper understanding of CCQE and CCl7r+ is essential for the reconstruction of the neutrino 

energy. 

The main task in a Vg appearance experiment like MiniBooNE is to detect electron neutrinos in a (almost) pure 
^ ■ beam. The signal event, the Vg CCQE interaction, is dominated by background. A major problem comes from 

misidentified events, mainly because of the fact, that the MiniBooNE detector cannot distinguish between a photon 
and an electron. Thus, induced neutral current tt" production, where the nP decays into two 7s, is the major source 
' of background when one of the photons is not seen or both Cherenkov rings overlap. 

As all of the present oscillations experiments use nuclear targets, it is mandatory to consider FSI, i.e., pion 
rescattering, with and without charge exchange, and absorption in the nuclear medium. A realistic treatment of the 
FSI can be achieved in the framework of a coupled-channel transport theory — the GiBUU model. 

After a brief review of our model, we first discuss the impact of pion production on CCQE measurements. We further 
investigate the influence of nuclear effects on CC 1 7r+ and NC 1 %^ cross sections, and, where possible, we confront our 
model to recent data measured at MiniBooNE. 
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GIBUU MODEL 

We treat neutrino-nucleus scattering as a two-step process. In the initial-state step, the neutrinos interact with nucleons 
embedded in the nuclear medium. In the final-state step, the outgoing particles of the initial reaction are propagated 
through the nucleus using a hadronic transport approach. 

In the energy region relevant for MiniBooNE, SciBooNE and K2K, the elementary vA^ reaction is dominated by two 
processes: quasielastic scattering and the excitation of the A resonance (P33(1232)). Additionally, our model includes 



the excitation of 12 A'^* and A resonances with invariant masses less than 2 GeV and also a non-resonant single-pion 
background. The nucleon, — A and N — N* vector form factors are based on recent electron-scattering data while 
the axial couplings are obtained assuming PC AC (partial conservation of the axial current). The dependence of the 
nucleon and N — A axial form factors, and Cj , is fitted to bubble chamber neutrino-scattering data. 

The neutrino-nucleon cross sections are modified in the nuclear medium. Bound nucleons are treated within a local 
Thomas-Fermi approximation. They are exposed to a mean-field potential which depends on density and momentum. 
We account for this by evaluating the above cross sections with full in-medium kinematics, i.e., hadronic tensor and 
phase-space factors are evaluated with in-medium four-vectors. We also take Pauli blocking and collisional broadening 
of the outgoing hadrons into account. Our model for neutrino-(bound)nucleon scattering is described in detail in 
Ref. [1]. 

After the initial neutrino-nucleon interaction, the produced particles propagate out of the nucleus. During propa- 
gation they undergo FSl which are simulated with the coupled-channel semi-classical GiBUU transport model (for 
details, see Ref. [2] and references therein). It is based on the BUU equation which describes the space-time evolu- 
tion of a many-particle system in a mean-field potential including a collision term. Nucleons and resonances acquire 
medium-modified spectral functions and are propagated off-shell. Herby we ensure, that vacuum spectral functions are 
recovered after leaving the nucleus. The collision term accounts for changes (gain and loss) in the phase-space density 
due to elastic and inelastic collisions between the particles, and also to particle decays into other hadrons. Baryon- 
meson two-body interactions (e.g., TtN — > ttAO are described by resonance contributions and a small non-resonant 
background term; baryon-baryon cross sections (e.g., A^A^ NN, RN NN, RN R'N, NN nNN) are either 
fitted to data or calculated, e.g., in pion exchange models. The three-body channels ttA'A' NN and AA'A' NNN are 
also included. The BUU equations for all particle species are thus coupled through the collision term and also through 
the potentials. Such a coupled-channel treatment is required to account for side feeding into different channels. 

CCQE AND CCl;r+ ENTANGLEMENT 

One challenge is to identify true CCQE events in the detector, i.e., muons originating from an initial QE process. To be 
more precise, true CCQE corresponds to the inclusive CCQE cross section including all medium effects, or, in other 
words, the CCQE cross section before FSL The difficulty comes from the fact that the true CCQE events are masked 
by FSl in a detector built from nuclei. The FSl lead to misidentified events, e.g., an initial A whose decay pion is 
absorbed or which undergoes "pion-less decay" can count as CCQE event — we call this type of background events 
"fake CCQE" events. We denote every event which looks like a CCQE event by "CCQE-like". 

In general, at Cherenkov detectors such as MiniBooNE CCQE-like events are all those where no pion is detected 
while in tracking detectors such as K2K-SciBar/SciFi CCQE-like events are those where a single proton track is visible 
and at the same time no pions are detected. 

To investigate the relation between the CCQE-like and true CCQE cross section, we show their ratio as a function 
of proton and pion momentum thresholds in Fig. 1. As the proton is not at all relevant for the CCQE identification in 
Cherenkov detectors, the ratio is independent of the proton momentum detection threshold (dashed line in left panel). 
This is very different in tracking detectors which rely on the detected proton — here the efficiency is reduced to « 10% 
at a proton momentum threshold of 0.5 GeV (solid line in left panel). Even at |/5 = the efficiency does not exceed 
80% because of charge-exchange processes that lead to the emission of undetected neutrons and because of secondary 
proton knockout that leads to multiple-proton tracks. Focussing on the right panel of Fig. 1 we find that the CCQE-Uke 
cross section increases for both detector types as |/?lthres increases. In this case even more events with pions in the final 
state appear as CCQE-like because these pions are below threshold and thus not detected. 

Fixing the flux normahzation with HARP's pion-production data, the MiniBooNE collaboration has presented their 
first, preUminary absolutely normalized total, differential and double differential cross sections for CCQE at this 
conference and finds an excess of about 35% compared to the total cross section measured at NOMAD, ANL and 
BNL [3]. We strongly emphasize that these absolute cross sections depend directly on the pion background subtraction 
which again is based on the Monte Carlo prediction (cf., Ref. [3]). 

In Fig. 2, we show our prediction for the double differential cross section at MiniBooNE in muon observables, all 
calculated with = 1 GeV. The left panel shows the true CCQE events. To emphasize the role of "fake" CCQE 
events, we show the ratio CCQE-like/true CCQE in the right panel. Unlike for monochromatic beams, the QE and A 
peaks are not distinguishable any more but strongly overlap. This fact makes a model-independent cut based on muon 
variables to subtract the background impossible. 
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FIGURE 1. Ratio of the CCQE-like to the true CCQE cross section as a function of the proton (pion) momentum detection 
threshold for CC on ^^C at Ey = I GeV. The solid lines are obtained using the tracking detector identification, while the dashed 
lines are for Cherenkov detectors. 
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FIGURE 2. Double differential cross section on '^C averaged over the MiniBooNE flux as a function of the muon kinetic energy 
and the muon scattering angle. The left panel shows the true CCQE cross section, the right panel the ratio of the CCQE-like to the 
true CCQE cross section. 



MINIBOONE'S CC1K+ MEASUREMENT 



In the left panel of Fig. 3 we give our results for the single-7r+/QE ratio for CC interactions on mineral oil CH2. The 
solid lines denote the CC 1 7i;+-like/CCQE-like result, the dashed lines stand for the true CC 1 TT+Ztrue CCQE result, and 
the dash-dotted lines give the vacuum expectation. Note that we have applied the Cherenkov detector identification 
criteria. 

We emphasize that nuclear corrections cancel out in the ratio, only as long as FSI are not considered ("true" 
vs. "free"). In general, the complexity of FSI prevent such cancellations as one can infer from the result denoted 
with "like" which does not coincide with the "true" and "free" ones. 

We further compare to very recent MiniBooNE data [4] (Fig. 3 left). Let us first focus on the data denoted with 
the triangles (upper data set). These are corrected for FSI using a specific Monte Carlo generator, i.e., they give the 
cross sections for bound nucleons "before FSI". As this procedure introduces a model dependence in the data, a fully 
consistent comparison is not possible. Ignoring this inconsistency, our calculation denoted with "true" should be the 
one to compare with. The agreement is perfect for energies up to 1.5 GeV, and still within their error bars for higher 
Ey The MiniBooNE data denoted with bullets (lower data set) is their result for the ratio of CCl7r+-like to CCQE-like. 
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FIGURE 3. Left panel: single- ;r+/QE cross section ratio for CC interactions vs. neutrino energy on CH2 together with recent 
data from MiniBooNE [4] (upper data set: corrected for FSI, lower data set: uncorrected for FSI). The solid lines denote the 
CCl;r+-like/CCQE-like result (Cherenkov detector definitions), the dashed lines stand for the true CCl;r+/true CCQE result, and 
the dash-dotted lines give the vacuum expectation, i.e., the sum of the nucleon cross sections (with two additional protons in the 
MiniBooNE case). Right panel: same but with a dipole form for Cj (see text). 



As these data are not corrected for FSI within a specific Monte Carlo event generator, this observable is less model 
dependent. Still, the energy reconstruction requires specific assumptions as well as the detector simulation. We find 
that our calculation clearly underestimates the uncorrected data. However, the perfect agreement with their corrected 
distribution indicates a significant difference between the pion absorption models. 

The underestimate of the pion/quasielastic ratio in particular at higher energies could be due to, among other 
possibilities, an underestimate of the pion production cross section or an overestimate of the CCQE-like cross section. 
Both, which means for the latter the fake CCQE events, depend directly on the input at the nucleon level, i.e., in 
particular on the axial form factor C5 of the A resonance (see Ref. [1] for details). In the above calculations, we have 
used 

[ci{Q^)\^^Ci{Q) 



1 , 




-2 

; (1) 



with a = -0.25, b = 0.04 GeV^, = 0.95 GeV obtained from a fit to the ANL bubble chamber data and the PCAC 
value of C5 (0) — 1.17. However, the ANL and BNL results differ with the BNL result being approximately 30% 
higher To account for the BNL findings, we introduce — following Ref. [5] — a dipole form factor 

[C^(e')],, = C^(0)(l+^) , (2) 

with M^^ = 0.94 GeV. The results for both parametrizations are shown in Fig. 4 and Fig. 5 together with the data of 
ANL and BNL. The shape of the BNL distribution is described well by both form factors (right panel of Fig. 4) while 
the ANL data are overestimated by the dipole form. Still one should bear in mind that there is a small non-resonant 
background [6] that has been neglected in this channel. For a discussion of possible Deuterium effects, in particular at 
lowg2^ seeRefs. [5, 7]. 

In the right panel of Fig. 3 we present our results for the single-7r+/QE ratio using the dipole form factor. We find 
a slight enhancement of the ratio but still underestimate the data. Note that the enhancement is much more moderate 
than in Fig. 5. This is due to the fact, that both, numerator and denominator, i.e., also the CCQE-like cross section, are 
increased when the dipole form factor is used and, as a consequence, the larger CCQE-like cross section compensates 
the enhancement in the pion cross section. From this comparison we conclude that an increase of the total pion 
production cross section on the nucleon compatible with the BNL data seems to be insufficient to describe this ratio at 
all energies. A similar result for the ratio has been recently obtained by Athar et al. [11]. 




FIGURE 4. Differential cross section Aa /dQ^ averaged over the ANL flux (left panel) and the BNL flux (right panel) [invariant 
mass cut at IV < 1 .4 GeV]. The data are taken from Refs. [8, 9]. The form factors are detailed in the text. 




Ev [GeV] 



FIGURE 5. Total CC pion production cross sections compared to the pion production data of ANL [Refs. [10] (•), [8] (■)] and 

BNL [[9] (x)] The solid line has been obtained with [C^(Q^)\, while [C^(e^)]„ was used for the dashed line. 



NCl;rO 

In the left panel of Fig. 6, we show our results for NC single-;r*' production off '^C as a function of the pion kinetic 
energy. We have averaged over the MiniBooNE energy flux which peaks at about 0.7 GeV neutrino energy [12]. In 
NC reactions the total pion yield is dominated by production, while ;r+ dominate in CC processes (for details, 
see Ref. [13] and references therein). Comparing the dashed with the solid line (results without FSI and spectral 
function vs. full calculation), one finds a considerably change. The shape is caused by the energy dependence of the 
pion absorption and rescattering cross sections. Pions are mainly absorbed via the A resonance, i.e, through nN A 
followed by IsN NN. This explains the reduction in the region around Tji = Q. \ — 0.3 GeV. Pion elastic scattering 
nN nN reshuffles the pions to lower momenta and leads also to charge exchange scattering into the charged pion 
channels. The vast majority of the pions comes from initial A excitation (dash-dotted line), their production in the 
rescattering of nucleons is not significant at these energies. 

The right panel of Fig. 6 shows the results for NC single- tt" production off '''O averaged over the K2K energy flux 
which peaks at about 1.2 GeV neutrino energy [14]. Compared to the left panel, the spectrum is broader and extends 
to larger T„ due to the higher neutrino energy. The reduction in the region around = 0. 1 — 0.3 GeV is mainly caused 
by the pion absorption via the A resonance (compare dashed and solid lines). Again, pion production through initial 
QE scattering is not sizable. 




FIGURE 6. Left panel: NC induced single-;!*' production on '^C as a function of the pion kinetic energy averaged over the 
MiniBooNE flux. The results were obtained with [C^(2^)]^. Right panel: same on '^O averaged over the K2K flux. The dashed 
lines show our calculation without FSI or spectral functions, both included in the full calculation denoted with the solid lines. The 
dash-dotted lines indicate the contribution from the A resonance to the full calculation. 



FINAL REMARK 



Overall, the impact of FSI on observables is dramatic and not at all negligible. A qualitatively and quantitatively correct 
treatment is thus of great importance. More model-independent data are certainly needed. 
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